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A monopole creation operator is constructed: its vacuum expectation value is 
an order parameter for dual superconductivity in that, if different from zero it sig- 
nals spontaneous breaking of the U(l) symmetry corresponding to monopole charge 
conservation. The operator is tested on compact U(l) gauge theory on lattice. For 
SU(2) gauge theory it clearly demonstrates that confinement is produced by dual 
superconductivity. 
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I. INTRODUCTION 

A possible mechanism of colour confinement in Quantum Chromo-Dynamics (QCD) is 
dual superconductivity of the vacuum |I]-f3]]. According to this scenario, the chromoelectric 
field is channelled into Abrisokov flux tubes ||, in the same way as the ordinary magnetic 
field is in superconductors: the word dual indicates the interchange of roles between electric 
and magnetic fields and charges. The chromoelectric field mediating the force between 
coloured particles is squeezed by Meissner effect into flux tubes of constant energy per unit 
length, giving rise to the confining linear potential. These flux tubes behave as strings . 
The existence of strings in hadronic physics is supported by phenomenology J7|||. They 
have also been visualized by numerical simulations of QCD on the lattice P,PH|. Some 
evidence in favour of dual superconductivity of the vacuum has been produced by Montecarlo 
simulations [|TT|. 



A clear cut test of the mechanism would be the detection of monopole condensation 
in the ground state, analogous to the condensation of Cooper pairs in the ground state of 
an ordinary superconductor. Condensation implies that the vacuum is a superposition of 
states with different charge, which, in turn, is nothing but a spontaneous breaking of the 
£7(1) symmetry related to charge conservation ||12|| . Such a breaking is signaled by a non- 
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vanishing vacuum expectation value (vev ) of any operator carrying non-trivial charge, as, 
e.g. the scalar field of the Landau- Ginzburg model of superconductivity. That vev is called 
a disorder parameter in the language of statistical mechanics. 

In QCD the monopoles which are expected to condense and generate superconductivity 
are Dirac monopoles of a residual {7(1) symmetry, which survives after a suitable gauge 
fixing, known as abelian projection An abelian projection is defined as the gauge trans- 
formation which diagonalizes any operator transforming in the adjoint representation of 
the gauge group. The points where two eigenvalues of that operator coincide correspond 
to singularities of the gauge transformation and have the topology of world-lines of point 
monopoles p|,p!5|. Such monopoles have been observed on the lattice Of course the 

location and the number of monopoles do depend on the choice of the operator used to define 
the abelian projection: a possibility [|IJ is that physics is independent of it. The relevant 
abelian degrees of freedom can also be fixed by a somewhat different procedure, known as 
maximal abelian projection |T5[ . Our operator allows to investigate unambigously what 
abelian projection, if any, defines the monopoles relevant to confinement. 

In this paper, we will present the construction of the operator, which is rather general 
and can be used for any kind of solitons. We will analyze how the operator works in compact 
i7(l) gauge theory |T7 . 



For that theory, there exists a construction of a disorder variable describing monopole 



condensation fL8 |, which is rigorous but based on a particular form of the action (the Villain 
action). Our operator coincides with the above one in the Villain case, but can be used with 
different forms of the action. This is particularly important in order to use it for non-abelian 
theories, where the effective {7(1) action after abelian projection is not known. 

For {7(1) we get a spectacular signal of monopole condensation [|l7j]. We then apply the 
same construction to SU (2) gauge theory: we consider the monopoles corresponding to the 
abelian projection in which the Polyakov loop is diagonal, and we investigate the condensa- 
tion of monopoles across the deconfining phase transition. Here again we find a spectacular 
signal of monopole condensation. The monopoles corresponding to that abelian projection 
do condense in the confined phase: Confinement is produced by dual superconductivity. 

The next steps of our investigation, which is in progress, will be to study the behaviour of 
monopoles defined by different abelian projections, to test t'Hooft ideas, and, in particular, 
the properties of the maximal abelian projection. 

In sect. 2 the construction of the monopole creation operator is presented. Sect. 3 contains 
the results for compact {7(1). Sect. 4 the results for monopole defined by abelian projection 
diagonalising the Polyakov loop in SU(2) gauge theory. Sect. 5 contains a few concluding 
remarks. 



II. THE MONOPOLE CREATION OPERATOR. 

For the sake of definiteness, we shall consider here only the case of a {7(1) monopole. 
Our procedure can easily be generalized to all kind of solitons. 

Let £? M (x, y) = (0, b(x, y)) be the classical field produced in the location x by a Dirac 
monopole at rest in y. We can make any choice for the gauge, e.g. by putting the string 
along the positive z-axis. Then, defining r = x — y and r = ||r||, we have: 
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(1) 



where g is the charge of the monopole, satisfying the Dirac quantization condition e g 
If Ilj(x, t) is the conjugate momentum to Aj(x, t), then the operator 



/i(y,t) = exp{i J d 3 x 6j(x,y) FL^x,*)} 



(2) 



creates a monopole at the location y and time t. This can be immediately seen in the 
Schrodinger representation of the fields, where we have: 



/i(y,0)|A(x,0)) = |A(x,0) + b(x,y)) 



(3) 



Equation (|3J) is a trivial consequence of the canonical commutation relations between the 
fields and their conjugate momenta, and is nothing but the field-theoretic equivalent of the 
familiar statement that e ipa translates the coordinate q by a. /i(y,t), applied to any field 
configuration, adds a monopole to it. This can be restated in terms of commutation relations 
as: 



[Ai(x, t), /i(y, t)} = bi(x, y)/i(y, t) 
[n,(x,t), / i(y,t)] =0 



(4) 
(5) 



which also show that the electric field of the configuration is left unchanged. It is worthwhile 
to notice that the specific choice of the gauge for b is irrelevant: what really matters here 
is topology, which is independent of it. Our operator [i is similar to operators introduced in 



the literature by different constructions, in various contexts [19|. 

Now, if the ground state of the theory has a definite monopole number N, then, under a 
magnetic 17(1) rotation, we have: U\0) = e ilpN \0). Since UfiW = e ilp fi, then (0|/x|0) = 
(0|/i|0)e^ which implies (0|/x|0) = 0. Therefore if (0\^\0) ^ then |0) is not U(l) invariant 
and there is spontaneous symmetry breaking of U(l). A translation by a static field g(x) 
such that curl g = in this language corresponds to a pure gauge transformation: 



7(t) = exp ji J <i 3 x gr*(x) IT(x, t) 



(6) 



By gauge-invariance (0|7(t)|0) = 1. Performing the Wick rotation to Euclidean space, we 
obtain: 



^e{y,Va) = expj-# j d 3 x 6i(x,y) IL,(x,?/4)J 



(7) 



From now on, we will be interested only in the Euclidean quantity and drop for simplicity 
the E subscript. We will first compute (/jl) = (0|/i£(y, 2/4) |0) for free photons. Rescaling the 
fields by a factor l/y/P with (3 = 1/e 2 , we have 



(H) = VA jexp -|y d 4 xF, v F, 
-(3 Jd 3 xF 0i (x, yi)bi(x, y 



(8) 



x exp 
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The integral is Gaussian and can be directly computed giving: 

(p) = exp 1 1 / ^<iV^(-A0>M-k) 6,-(k)J 

(^(fc)F 0j (-fc)) = ^,- k2 ^~ k ^ (9) 
If b is such that kibi(k) = 0, we have, performing the integral over ko in the second term, 

<">=« p f {/(i^ |b(fe)|2 -|/(0^ |k||b(fc)|2 } (10) 

In the analogous calculation for a gauge transformation 7, <7i(k) oc fcj, the second term of 
eq.(|S]) does not contribute, and 

w=^{i/S |g( * )|a } 

We realize that the first term in the exponent of eq.(|I]) is a normalization which can be 
subtracted by taking instead of the ratio (p) = (p) / (7) where (7) is defined by means 
of any gauge transformation g, such that: 

' (*or= 44-Mk)\ 2 (12) 



(2tt) 4IOV n J (2tt) 4 1 

Then for free photons, which, in the compactified version of the theory, correspond to the 
deconfined phase, j3 ^> (3 C 

^)=ex P {4/7^|k||b(fc)| 2 ) (13) 



A J (2tt) 31 

The integral in the exponent, once regularized at small distances, tends to +00 as V — > 00. 
In the infinite volume limit (p) = as it should be, since the perturbative vacuum has 
zero magnetic charge. The same situation appears when computing the overlap of the Fock 
vacuum to the Bogolubov rotated vacuum in a superconductor PD1 . 

An alternative way of looking at the normalization factor in eq. fllCf ) is to go back to the 
very definition of Feynman path integral: when computing the vacuum expectation value of 
an operator like e vp ^ a , after the usual discretization of time is performed in intervals of size 
S, the operator appears in a matrix element of the form 

iHS ipa\ v f dp . 1 iHS ipa 



(x n+1 \e-* Hd e tpa \x n }= / -J-(x n+l \e-* Hd e^ a \p)(p\\x n ) 

(2tt) 



dp - 



(27T) 

The integral over p can be performed giving 

e -i(x n+1 -x n ) 2 ^+ia(x n+1 -x n )m+V (x„) ^-ia? 

The first factor is the lagrangian definition, the second factor corresponds to the subtraction 
operated in going from (p) to (p) in eq. (|1"3|) . 
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III. LATTICE FORMULATION. U(l) GAUGE THEORY. 



A lattice version of the operator p is obtained by replacing eF 0i by the plaquette P 0i , or 
better, by its imaginary part: eF 0i — > ImP 0i and discretizing the field 6j. Then the disorder 
parameter becomes: 

(li) = ^fvUexp{-P[S+ £ b t (n,y)ImP 0l (n)]} 

or, if we want to cancel the unwanted normalization, we can divide by: 

(-Y) = ^JvUexp{-/3[S+ £ gi (n)ImP 0l (n)]} 

obtaining 



no=yi 



<fi) 



JvUexpl-P[S+ ^6 8 (n,y)/ m ? 0! ( n )] 



n=y 4 



JvUexpl-f3[S + Y,9i{n)ImP u {n)\ 



(14) 



n=y 4 



We stress once again that any gauge function is acceptable, provided the normalization 
condition ( |T2D is satisfied. In the following we shall use Wilson's action S = ^ Ym,u,v(^~Pia>)- 
If we blindly compute (p) or (p) by numerical simulations, a first technical difficulty 
arises. We are faced with the usual problems encountered in computing quantities like 
a partition function, which are exponentials of extensive quantities, proportional to the 
number of degrees of freedom. The distribution of the values is not Gaussian and the error 



does not decrease by increasing statistics (see, e.g. pi |, where the same problem appears in 
a different context). To avoid that, we will compute the quantity: 



^log(/,) = ^log(/,)-^log( 7 ) 



(15) 



At P = 0, (/i) = (7) = 1, and therefore: 



(jj) = exp 



dP'p(P') 



Putting S b = Y,bi{n)ImP 0i (n) 



"o=0 



and S g = J29i( n ) ImPi 



0i lnn=0 



we get: 



P — (S + Sg) S + Sg - (S + S b ) S + S b 



(16) 



(17) 



which can be evaluated by numerical simulations. The subscript on the average indicates the 
action defining the Feynman integral. The two quantities on the rhs have the same strong 
coupling expansion. Thus, the use of (fl) instead of {p), besides producing a cancellation of 
the spurious normalization of the Feynman path integral, can help in eliminating the lattice 
artefacts produced by the discretization which can spoil the continuum limit. This brings us 
to a second, more physical difficulty, which is the continuum limit. In fact, while for QCD, 
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which is asymptotically free, we expect that, at sufficiently high (3, lattice artefacts should 
cancel, for a model like £7(1) this point is not so clear in principle. In Ref. a proof is 
given of monopole condensation in the confined phase of U(l), defining a disorder variable 
for the Villain action. We have checked that our (p) operator exactly coincides with the one 
of Ref. |18j, when we use the Villain action: we expect that for Wilson action the same will 
hold. We have then computed numerically p. 

For a good order parameter, we would expect p to be zero, or (p) = 1 below the critical 
value of the coupling (3 C and then to show a large negative peak around f3 c , corresponding to 
a drop to zero of (p). At larger values of [3, we have free photons and Eq. (O) should hold. 



Figure 1 shows the behaviour of p for a 12 4 lattice, for a monopole in the center of the 
space lattice. In order to be able to identify the signal as a genuine physical result (i.e. not 
due to lattice artefacts), we have performed a number of checks: 

1 . We have changed the form of hi by a gauge transformation to get the Wu-Yang expression 
of the monopole potential. The result does not change qualitatively. 

2. (7) shows practically no signal at f3 c within the errors, and does not change appreciably 
by changes of gi. To test that we have computed (7) for two different choices of g,(x), both 
satisfying the condition Eq. (|T2"|) : (71) for gi(x) = const, and (72) for g 2 (x) oc x/|x| 2 . In 
Fig. 2 we display p gauge = ^ In |*4 and p = ^ In -^K for a 6 4 lattice. p ga uge shows no relevant 
signal at (3 C . 

3. We have measured the correlation function between a monopole antimonopole pair 
at large time distance, in the same position in space. To do that we define C(d) = 
(p(0,0) p(0,d)) and 

S b M = £ 6,(11, 0) [ImP oi (n, 0) - ^(n, d)} 

n 

We measure ^ In ^ = p (b ; b) (d) or 

P(b,b)(d) = 2(S + S g )s+s g - (S + S bb -(d)) s+ s bi - (S) s 
1 we have 



Since ^ 



(3=0 



C(d) 



rp 

J o P(b,b)(d)d{3 



(18) 



By the cluster property we should have at large d that C(d) — > (p(0, 0)) 2 , or p^ t b)(d) — > 2p. 
Fig. 3 shows that this expectation is indeed verified. We notice that the height of the negative 
peak of p at f3 c increases with volume [see Fig.l and Fig. 2]. The value of f3 c as defined by 
the position of our peak is j3 c = 1.01(1) for a 6 4 lattice and j3 c = 1.009(1) for a 12 4 lattice. 

We have taken monopole charge n = 4 to get a good signal with a relatively low statistics 
(tipically 10 4 configurations per value of (3): smaller charges (n = 1,2) give similar results 
but the signals are smaller and more noisy. Our statistical errors are shown in the figures, 
when they are larger than the symbols used. 
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IV. SU(2) GAUGE THEORY. MONOPOLE CONDENSATION AND 

CONFINEMENT 



The monopoles which should condense in QCD vacuum in the confined phase are Dirac 
monopoles of the U(l) field defined by the abelian projection The density of current of 
such monopoles is 



3u = d^F t 



where, in the notation of ref. [22 



F, v = Xo a G% - ^-e ahc Q a {D,Q b ){D u Q c ) (19) 

IQI g\Q\ 

a,b,c— 1, 2, 3 

Q a is the Higgs field in the Georgi - Glashow model. Abelian projection is a gauge trans- 
formation which brings Q in the third direction, Q = (0, 0, Q 3 ). 

In QCD the role of Q a should be played by some operator trasforming in the adjoint 
representation |IJ: what this operator is, if any, is the problem under investigation. 

After abelian projection the second term of eq.(^) gives zero and F^ v = G 3 . 

As a first candidate for Q a we have chosen p a = |Tr(— i In P a a ) where P is the Polyakov 
line exp(J A o dx ). This choice corresponds to take Aq as the Higgs field. For that choice 
D Q = and 

v a 

Since the expression (|19D is gauge invariant, we do not need to perform any gauge transfor- 
mation to make InP diagonal. We simply define the gauge invariant quantity 

The operator creating the relevant monopole is 

At(y, t) = exp{- J rf 3 x6i(x, y)F M (x, t)} (20) 

Exactly as for U(l) we can define (p) = (fj) / (7) where 

( 7 (0 ) = exp{-| rf 3 x^(x)F 0i (x,t)} 

We shall take for ^(x) a constant field, such that eq. (|T2|) is satisfied. 
The lattice version of Fq{ will be 

l Tr{-2 lnPn ,} 
01 2 \p\ 

and 
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fouexpl-f3[S+ $>(n,y)F 0j (n)]} 

(A> = ~ ( ^ f (21) 

] T?f/exp + £ ^(n)F 0i (n)] 



n=j/ 4 

We will define p as in eq.flTSD, and therefore again: 



(fi) = exp 



o 



d0p{0 



As in U(l) we measure the behaviour of p{(3). across the deconfining phase transition The 
simulations are performed on asymmetric lattices, with temporal size smaller then spacial 
size N s . 

Fig. 4,5,6,7 show the behaviour of p for 8 3 x 4, 12 3 x 4, 12 3 x 6 and 16 3 x 6 lattices. The 
behaviour of p is similar to the U(l) case: p is compatible with zero below the deconfining 
transition, has a sharp negative peak at the deconfining phase transition, and a negativetail 
at high p, corresponding to a phase of free "photons" (see eq.flT^p. Qualitatively the shape 
of (p) is 

(p) = m - p) 

The negative peak increases with space volume at fixed Nt (compare 8 3 x 4 with 12 3 x 4, 
12 3 x 6 with 16 3 x 6) which means that as V — > oo, for (3 > /3 C , (p) — >■ 0. 



The peak is placed exactly at the official location of the phase transition for Nt = 4 p3 
and moves as expected from renormalization group from Nt = 4 to Nt = 6. 

The usual order parameter | (P) | is plotted for comparison in the figures. 

The depth of the peak does not change appreciably from 8 3 x 4 to 12 3 x 6, showing that 
the signal mainly depends on the physical volume of the lattice, i.e. on its size in fm 3 , and 
not on the geometrical volume of it. 

We have also measured the correlation function 

(/2(10,0,0,0),/2(0, 0,0,0)) 

on a 8 2 x 20 x 4 lattice to check the cluster property, and tested that it is equal to the square 
of the disorder parameter on the same lattice (Fig. 8). 

We can finally conclude that the monopoles defined by abelian projection identified by the 
Polyakov loop are charges of an £7(1) symmetry which is spontaneously broken in the confined 
phase, and restored above /3 C , or that confinement is produced by dual superconductivity. 



V. CONCLUDING REMARKS 

We have demonstrated that confinement in SU(2) gauge theory is due to dual super- 
conductivity of the vacuum. We have identified the monopole charges which condense and 
produce it. Our method consists in a direct detection of spontaneously breaking of the U(l) 
symmetry corresponding to monopole charge. 

The method has been successfully checked on U(l) compact gauge theory. 

A number of open questions are currently under investigation 
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Are there different abelian projections defining monopoles which condense in QCD 
vacuum? Is t'Hooft's idea about the physical equivalence of different abelian projec- 
tions correct? In particular we are investigating the so called maximal abelian gauge, 
defined by the gauge transformation which maximises the quantity 

J2 Tr {U^(n)a 3 UjXn)a 3 } 



On this gauge there is a lot of evidence collected in favour of dual superconductivity 
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The extension to SU(3) of the above construction. 

The application of our method to other systems, like e.g. the 3d x-y model, in which 
a phase transition is driven by the condensation of solitons. 
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FIGURE CAPTIONS 



Fig.l p versus (3 on a 12 4 lattice 

Fig. 2 p and p gaU ge versus (5 on a 6 4 lattice. 

Fig. 3 Pbi(d) versus /3 at d = 4, 7, 9, compared to 2p. 

Fig. 4 p (circles) and (\P\) (squares), 8 3 x 4 SU{2) gauge theory. 

Fig. 5 p (diamonds) and (\P\) (squares), 12 3 x 4 SU(2) gauge theory. 

Fig. 6 p (triangles) and (\P\) (squares), 12 3 x 6 SU(2) gauge theory. 

Fig. 7 p (open squares) and (\P\) (squares), 16 3 x 6 SU(2) gauge theory. 

Fig. 8 p b i(d) at d — 10 (circles) compared to 2p (triangles) and their difference (squares). 
SU (2) gauge theory. 
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